www and Internet models from 1955 till our days and 
the "popularity is attractive" principle 
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We note that the model discussed in the commu nication of S. Bornh oldt and H. Ebel (World Wide 
Web scaling exponent from Simon's 1955 model, 



model considered and solved exactly in our paper, cond-mat/0004434. These models may be used 



ond-mat/0008465) is the particular case of the 



for estimation of the order of the deviation of the scaling exponent from 2 both for the distributions 
of incoming links and links coming out from nodes but not for the obtaining some specific values of 
the exponents from the WWW growth data. We emphasize that, unlike the statement of Bornholdt 
and Ebel, both the network under consideration and the model of Barabasi and Albert provide quite 
equal possibilities for individual growth. There is no great difference between them in this respect. 
The resulting distributions for in dividual nodes and arising scaling relations have been obtained 
in our paper, :ond-mat/0004434. We discuss briefly the modern state of art in the physics of the 
evolving networks and the great role of the general principle - popularity is attractive - in the self- 
organization of complex communications networks, in physics of nonequilibrium phenomena, and in 
Nature. 



An interesting linking between the Simon's model ||] 
(1955) and new conceptions of the developing networks 
was made in the communication of Bornholdt and Ebel 
. Using a simple model of the growing network with 
preferential linking and referring to the Simon's model 
Bornholdt and Ebel have obtained the connectivity dis- 
tribution P{k) oc k~'^ with the exponent 7 which value 
have coincided with the experimental result (2.1) [g[Q] 
for the exponent of the distribution of incoming links in 
World-Wide Web. 

Here, we note that this model is the particular case 
of the model which was considered and solved exactly 
(for large networks) in our paper ||5|]. The only differ- 
ence is the different definitions of the time scales that 
do not influence the result. We show that these models 
may explain the scales of the deviation of the exponent 7 
from 2 for both the distribution of incoming links in the 
World Wide Web and the distribution of links coming 
out from nodes but they can not provide the particular 
experimental values of these exponents. We demonstrate 
these estimations and explain their quality. 

We point out that the statement of Bornholdt and 
Ebel that this model "allows for different growth rates" 
for nodes with same connectivity and, e.g., the model of 
Barabasi and Albert Q do not allow is not right. Both 
models provide quite equal possibilities in this respect. 
Their connectivity distributions of individual nodes are 
of a similar type. 

We describe the connections of the idea of the preferen- 
tial linking with classical problems of statistical physics of 
nonequilibrium phenomena and discuss briefly the state 
of art in the statistical physics of evolving networks. 

Also, the aims of the present comment are to explain 
popularly the results of our papers nQM obviously writ- 



ten in a too mathematical language to be noticed and to 
manifest quite obvious simple ideas. 

Below we: 
(i) describe the model and its particular cases, 
(ii) explain what estimates for the World Wide Web ex- 
ponents it may give, 

(iii) explain what kind of the distributions of connectivity 
of individual sites does the preferential linking produce, 
(iv) discuss the modern state of art in understanding of 
the World Wide Web, Internet, and other evolving net- 
works, and, finally, 

(v) demonstrate the simplest derivation of the expression 
for the exponent of the connectivity distribution, Eq. (Q). 

/. The model. — At each increment of time a new node 
is created. Simultaneously m directed links is distributed 
among all the nodes of the network by the following rules: 

(i) The source of any of these links may be anywhere; 

(ii) The target ends of these links are attached to nodes 
chosen preferentially: Probability to chose some particu- 
lar node is proportional to A: -I- ^. Here, k is the number 
of the distributed preferentially links which point to this 
node, and A is some positive constant. 

The distribution of incoming links in this model is 
of the power-law form with the exponent 7 — 2+A/m [p|. 
We call A the initial attractiveness of nodes (for links). 
A is a parameter of the problem and is any positive con- 
stant. In particular, if n links is already directed to any 
newborn node, A = n + B, where B is another constant. 

We stress that the preferential linking with the proba- 
bility proportional to fc -|- const is the general type of the 
preferential linking producing free-scale networks (i.e., 
the networks with power-law connectivity distributions). 
E.g., the rule with the linking probability proportional 



to k"', a ^ 1, do not give scale-free networks [Q. Hence, 
we consider the most general situation. 

In addition to links distributed preferentially, at each 
increment of time, one may also distribute some links, 
Tir, randomly without any preference (always there ex- 
ist some crazy guys which make their references to God 
knows what). Hence, 

— TO is the number of links distributed preferentially each 
unit of "time" , 

— n is the number of links pointed at a new node at the 
instant of its birth, n > 0, 

— rir is the number of links distributed without prefer- 
ence at each increment of time, rir > 0, 

— i? is any constant such that n + i? > 0, 

— one node is born per unit of "time" . 

With these definitions, the expression for the exponent 
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Note that all the factors accounted for in Eq. p) (see the 
numerator), in fact, play the same role. An exact form of 
the connectivity distribution was obtained in ^. Below 
we show how our relation, Eq. (nh, may be obtained in 
a most simple way using continuous approach. We have 
demonstrated that this approach gives exact results for 
the scaling exponents (5|j^,P] . 

Bornholdt and Ebcl considered the model in which, at 
each increment of time, one link is added to a network. 
In our case, time is equal to the total number of nodes. 
If we rescale time and set n^ = 0, n = 1, and B = 0, 
we get the particular case considered by Bornholdt and 
Ebel with their probability a = 1/(1 -|-to) (one may check 
that all the quantities may be taken noninteger). If we 
set rir = and n + B = m, we get, in fact, the Barabasi- 
Albert's model. 

//. Possible estimations for the World Wide Web. — We 
do not know values of any of the quantities from Eq. (Q). 
The constant B may takes any values, the number of the 
randomly distributed links, rir, may be not small (many 
crazy guys), n is not fixed. From the experimental data 
Q, we know more or less the sum m-|-n-|-nr~10»l 
(between 7 and 10, more precisely), and that is all. 

The only thing we can do, it is to fix the scales of 
the quantities. The natural characteristic values for 
Ur + n + B in Eq. (|I|) are (a) 0, (b) 1, (c) m > 1, 
and (d) infinity. In the first case, node have zero initial 
attractiveness, and all new links are directed to the old- 
est node, 7 — > 2. In the last case, there is no preferential 
linking, and the network is not scale-free, 7 — > 00. Let 
us consider the really important cases (b) and (c). 

(b) How do pages appear in the Web? Suppose, you 
want create your personal home page. Of course, first 
you prepare it, put references, etc. But that is only the 
first step. You have to make it accessiable in the Web, 
to lanch it. You come to your system administrator, he 
put a reference to it (one reference) in the page of your 



institution, and that is more or less all - your page is in 
the World Wide Web. If the process of appearing of each 
document in the Web is as simple as the creating of your 
page - only one reference to the new document (n = 1) - 
and if one forgets about the terms Ur and B in Eq. (Q), 
than, for the exponent of the distribution of the incom- 
ing links, we would get immediately the same estimation 
as in (2), 7 — 2 ~ 1/m ~ 10~^. This estimation is in- 
deed coincides with experimental value 7™ — 2 = 0.1 [Q. 
We repeat again, that it follows only from the fixation of 
scale of the involved quantities. We emphasize that there 
are no any general reasons to set, e.g., B = 0. A lot of 
real processes are not included in this estimation. Ag- 
ingof nodes changes 7 mKV , account for dying of nodes 
7 [H (the half-life of a page in the Web is of the order 
of half a year) changes 7. The ratio between the total 
number of links and the number of nodes in the Web is 
not constant [H, it increases with time, the growth of the 
Web is nonlinear. This factor also change the value of 7, 
in future it may become even lower than 2 [021. 

(c) Above we discussed the distribution of incoming 
links. Eq. (^ may be also applied for the distribution 
of links which come out from documents of the Web. In 
this case all the quantities in Eq. (R) takes other values 
which are again unknown. Nevertheless, one may think 
that the number of the links distributing without any 
preference, n^, is not small now. Even beginners proceed 
with linking of their pages. Hence, n^ ~ m - we have 
no another available scale, - and 7 — 2 ~ m/m ~ 1. 
We again compare this estimation with the experimental 
value, "font — 2 = 0.7 Q and find that we have made only 
small mistake. Of course, this is not quite fair but now 
we see what estimations may be made. 

What may be done to improve these estimates? As we 
saw, any experimental information about the correlation 
between the numbers of in- and out- links of nodes would 
be useful. 

///. Scaling of the distributions of connectivity of indi- 
vidual nodes. — One may study not only the global con- 
nectivity distribution of the network but the connectivity 
distributions of individual nodes. We have obtained the 
exact expressions for these distributions of the model un- 
der discussion Q . Unlike the statement of Bornholdt and 
Ebcl that this model provides some extra possibilities for 
different growth rate, we found that these distributions 
are of the same type as in the model of Barabasi and Al- 
bert H . It is a mistake to think that the network under 
consideration may demonstrate the growth crucially dif- 
ferent to the Barabasi- Albert's model. Indeed, the last 
one is only the particular case of the introduced model 
but this is not a marginal case to produce a radical dif- 
ference. Therefore, the "rich-get-richer" behavior also 
present in the model under consideration, although all 
these networks provide broad distribution of the follow- 
ing scaling form: 



p(fc,s,t)=(-j / 



1 



(2) 



Here, p(k, s, t) is the connectivity distribution of the node 
appeared at time s measured at time t (i.e., when the size 
of the network is i), s < t. The scahng exponent (3 is re- 
lated to 7 by the following scaling relation: 



/3(7-l) = l, 



(3) 



that was obtained in 0. It was proved in [gj, that this 
relation is general and is valid for all such scale-free net- 
works. We found that the scaling function is broad. 



/(^) 



"^exp(— x), so the connectivity distributions 



of individual nodes are broad. They are narrow only in 
the continuous approximation. 

IV. The principles of communications networks. — 
Perhaps, the greatest results in the theory of communi- 
cations networks belong to "the farther" of the Internet, 
P. Baran, (1964) pSl, so they were obtained long ago. 
Nevertheless, the idea of the preferential linking in appli- 
cation to the evolving networks is really strong. Indeed, 
one of the first questions concerning the communications 
networks - why are they scale-free? 

Two answers have been given quite recently. First, 
- because of the preferential linking |g]. The growing 
networks are self-organized into the scale-free structures 
because popularity is attractive. This general principle 
governs the process of the growth of the World Wide 
Web! 

Second, — because otherwise they would be unsta- 
ble, weak against processes of decay, and could not ex- 
ist as united systems mm. It has been shown that only 
the scale-free networks with the exponent 7 less than 3 
are resilient to random breakdowns |l£|. The "infinite" 
scale-free networks with 7 < 3 do not decay for any con- 
centration (less than one) of randomly removed nodes or 
links. 

Therefore, the evolution of the World Wide Web is 
governed by one of the most general principles in Nature. 
Only lazy can not imagine its numerous realizations in 
natural and social sciences, and in economics. It certainly 
establishes not only networks of movie actors. In partic- 
ular, in physics, it produces various processes of growth 
and fragmentation. For instance, in g we studied the 
growth of the networks mathematically as the problem 
of the preferential distribution of new particles among 
growing number of places. 

Now the World Wide Web is the most striking product 
of the principle popularity is attractive but it is only one 
of its products. It is wonderful that even in 1955 people 
thought about the quantitative formulation of this global 
principle. 

V. Derivation of our expression for the exponent. — Let 
us show in the simplest way how one may obtain Eq. M). 
We use the continuous approximation which gives exact 
results for the scaling exponents H . Using this approach, 
one may write the following equation with the bound- 
ary condition for the average connectivity (here, average 
number of incoming links) of the node s at time t: 



dk(s,t) Ur k{s,t) + B 

— 'rn 



dt t J*du[k{u,t) + By 

k{t, t)^n. (4) 

The first term in the right hand part of it arises from 
the links distributed without preference, the second one 
- from the preferential linking. Integrating Eq. (H) from 
to i we get 



* , dk(s,t) d f\ -, , -, , 
ds — jr-^ ^ Y I dsk{s,t)~k{t,t) ^Ur + m (5) 



and the obvious relation. 



ds fc(s, t) = [Ur + 7TT. + n)t . 



(6) 



Substituting Eq. (^ into Eq. (Q) we obtain the equation 
dk{s,t) Ur ni k(s,t) + B 



dt t m + Ur + n + B t 



(7) 



the solution of which in the scaling region is of the fol- 
lowing form: 



fc(s, t) oc 



S\-/3 



(8) 



(we used the boundary condition k{t, t) = n) with the 
scaling exponent 



P- 



m + Ur + n + B 



(9) 



Finally, using the scaling relation, Eq. (^, we obtain 

7=1 + 1//3 = 2+ {nr + n + B)/m. 
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